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Abstract 



The author tries to derive the asymptotic expression of the large eigevalues 
of some vectorial Sturm-Liouville differential equations. A precise description for 
C^ ■ the formula of the square root of the large eiegnvalues up to the 0(l/n)-term is 

obtained. 
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1 Introduction 

For the scalar case, it is well-known that: if q^^'{x) G £""^(0, vr), the equation 

-y" + q{x)y = Xy, y'{0) - hy{0) = 0, ^'(vr) + Hy{Tr) = 
has eigenvalues possessing the following asymptotic form 

^ " ^ n n2[™/2]+i ^ n2[W2]+i' 

where Uo = {h + H + 1/2 J^ q{t)dt)/'ir and the other Oj, 1 < j < [m/2] are 
constant depending on h, H, q{x) and the higher order derivatives of q{x) up to 
order m. The main technique for deriving above expression is by using iterative 
method. One can find the standard approaches in the appendix in Q or in |^]. 
For the vectorial case, the author only have seen analogous results in 0] for the 
case d = 2 and the boundary condtions are different to the cases the author has 
considered. Motivated by the arguments in [^, we still can use iterative method 
analogous to the one in [^, but only partial result as Theorem [l] can be obtained. 
The technique we use in this section is similar to the one in Q, but the major 
difference between theirs and the author's is that the author locate the eigenvalues 
by using the associated matrix differential equation and by the theory of complex 
operator-valued functions. Hence, we can dispose any dimensional case. Consider 
the following vectorial Sturm-Liouville differential systems 

-0"(x)+P(x)0(x)=A</.(x), cl)'{0) + HL(l>iO)=0 = <P'i7T)+HR<P{7r) (1) 

where P{x) is an N x N real symmetric matrix-valued function and Hi, Hr are 
N X N real symmetric constant matrices 

Consider the matrix differential equations related to (|l|) and (FS), respectively, 

-Y"{x) + P{x)Y{x) = XY{x), y(o) = /, y'(o) = -Hl, 
-C"{x) = XC{x), C(0) = /, C'(0) = -Hl, 

the solutions Y{x; X) and C{x; X) are connected by the following identity 



Y{x;X)=C{x;X)+ IC{x,t)C{t;X)dt, VA (2) 

Jo 

and C{x; A) = cos(vAx)I — sin(vAx)i:fL/vA. The connection (||) is easy to prove, 
one can see a more general version in [^ . 

In this paper, we also use the symbol {P, Hl, Hr) to denote the eigenvalue 
problem (||). Denote Ti{P,Hl,Hr) the sequence of eigenvalues of {P,Hl,Hr). 
We can arrange those elements in it in ascending order as 

/^o < /Wi < • • • < /^fc < • • • • 



We use the symbol cr{P, Hi, Hr) to represent the set of eigenvalues of (P, H^, Hr), 
and arrange those elements in this set in ascending order as 

Ao < Ai < . . . < Afc < . . . . 

The multiplicity of each eigenvalues A^ is denoted by m^. For each positive integer 
n, an index set A„ is defined as 

An = {k: A^/' G 5i/4(n)} 

where Br{£,) = {z : \z — 6.\ < r}. 

With above basic terminology and notations, we shall give the proof of the 
following theorem in section 2: 

Theorem 1 Suppose Xk G a{P, Hl^Hr) with multiplicity mk- For n is suffi- 
ciently large, if k £ A„, then 

A|=n+^ + 7n(fe), (3) 

where Ok is the characteristic value of the NxN real symmetric matrix ^(| Jq P{t)dt+ 
Hji — Hi), and 7n(fe) = o(l/n^). Furthermore, J2k€A„ "^fc = ^, where N is the 
dimension of the system (yj. 

2 Asymptotic analysis of eigenvalues 

In order to make this section become more readable, the author outlines the sketch 
of the whole approach as below: The eigenvalues A^ , /c > of (||) can be locate 
by determining whether the matrix-valued function 

W{Xk) = Y'{7r;Xk)+Hj^Y{7r;Xk) 

is singular or not. We use iterative method so as to write 

P^(A) = ^(A)+£:(A) (4) 

where we shall see that ^(A) has a quite neat and simple form from which we 
can determine those values A making ^(A) be singular. By the extension theorem 
of Rouche's theorem on operator-valued functions, we are getting close to locate 
the eigenvalues of (||), although we can not locate them completely and explicitly 
as we did in the scalar case. With above short sketch in mind, we now start the 
lengthy and tedious computations on the asymptotics of eigenvalues of (|l|) . 



Rewrite 



Y{x; A) = C{x- A) + T ^^^^^l^^Jlp{t)Y{t■, X)dt. (5) 

Jo vA 



Then 



Yi{x;X) = Y{x;X)-C{x;X) 

= Ii + f ""^""^^ Pm(t;A)dt. 

Using integration by part, we can compute /i as 

sinVAx,^, , cosvAx,,^, , ^^ 1/^/ n „/ nxn ^,,,, Sn 
h = J^ IC{x,x) + ^{}C{x,x)Hl + -{P{x) - P(0))) +0(|A|-2). 

Combining with (^), we have 

Fi(x;A) = '-^I^IC{x,x) + ^^^^{IC{x,x)HL + \{P{x)-Pm 

+ /; ""^""^^ Pm(t;A)dt + 0(|A|-|). 

As |9A| < K for some constant k fixed, by Gronwall's inequahty, we can prove 
that ||Yi(x; A)||oo = 0(|A|^^/^), as |A| is sufficiently large. From above argument, 
we have 

Y{x;X) = cos^/Xx{I + \lC{x,x)HL + -r{P{x)-P{0)) (7) 

A 4A 

On the other hand, from (^, 

Y'{x; A) = C'{x; A) + / cos VX{x - t)P{t)Y{t; X)dt. (8) 

Jo 

Plugging (|7[) into (g), and repeating the analysis as above, we have 
Y'{x;X) = —vXsinvXxI + cos vXx{}C{x,x) — Hl) 

sin vAx ,P(x) + P(0) 1 f^ -r^, .,^, ., ,^/ . ^^ s ,„s 
+ ^^( ^ \ +2i) P^*^'^(^^^)dt-}C{x,x)HL) (9) 

+ "-^^{\{p'ix) - p'm + \{[ pm{t,t)dt)HL 

+ \ J\p\t) - P{t)Pmdt + l{P{x)IC{x, X) - Hl)) + 0(|A|-i) 



By (^) and @, we have 



sm vAvr cos vAvr ^ ^ , , , _ 3 , 



T^(A) = -^/Asin^/A^/ + cosyA7rgi+ ^ ga + , g3 + 0(|Ar^), (10) 

V A A 



where 



gi=HR-HL + ICi7T,7T) 

P{7r)+P{0) , 1 



Pl-jr) + P(0) 1 /■'^ 
6^2 = ^ '^ +o / P(t)/C(t,i)dt + ^«/C(7r,7r) 

-}C{7T,7r)HL-HRHL, 



and 

p'fvrl - P'fO) 1 /■'^ 1 

^3 = ^^^ +2V0 P(t)IC{t,t)dt)HL + -Hn{P{n)-Pm 

+ ^ r(p2(t) - P{t)P{0))dt + ^(P(7r)/C(^, tt) - Hl) + ^ij/C(7r, ^)i?L 
o Jo ^ 

Write A = fp. Then the term \I'(A) in (^ is defined as ^(^^) = — //sin/i7r/ + 
cos^TT^i. The remainder terms of (|l^) is denoted by <?(/U^). Denote 'uj{^) = 
det(VF(;U^)) and fi = a + it. By the following identities 

sin(cr7r + iivr) = sincrvr cosh tvr + isinhtvr cos crvr, 
cos(o"7r + itvr) = cos crvr cosh tvr — i sinh tvr sin crvr 

and using the Laplace expansion of determinants, we have the following result: for 
any complex /x, vudj.) = 'i?(/x) + e{fi) where '!?(/u) = fi^ sin fnr and |e(/u)/'i?(/u)| = 
0(1/1^1). With above facts and Rouche's theorem, there lie exactly A^ zeroes of 
vadj.) within a suitable neighborhood of any sufficiently large integer n. If we 
denote a zero of Tu{^i) by n + 6ni, it is not difficult to see that 6ni = 0{l/n). In 
fact, we can compute lim„^oo ^ni, and hence we have the first part of Theorem 

To touch down the goal, let us back to (^). The matrix- valued function ^(/x) 
actually provides us much information. Since the matrix Qi is real symmetric, 
there exists an orthogonal matrix U, U~^ = U*, such that 

U*^{lj)U = diag[£ii(;u), Q2{^j), • • • , Qn{ij)] 

where ^i(/i) = — ^sin/ivr + ctj-ttcos^vt, 1 < i < A^, and o"j, 1 < i < iV are those 
characteristic values of constant matrix Qi. By (HI) and (p^, we have 



U*W{^i^)U = U*'^{fi^)U + U*£{fi'^)U (11) 

11 n^ 



u*^{^,')u + '-^^^J^u*g2U + o(^). 



Using the following Rouche's Theorem for operator- valued functions in [^, we 
can suitably locate the eigevalues of (Q) from (|ll|) . 



Theorem 2 Let W{fj,) and S{fi) be C{7i)-valued holomorphic function defined 
on a region 0, enclosed by a contour T. Suppose W is normal with respect to T 
and letV = W + S. If 

\\W-Hf,)S{fi)\\<l, (12) 

then V is also normal with respect to T, and 

m{r;V{-)) =m{r-W{-)). 

One can refer to [^, chapter XI, §9, 205-211] for the definition of all the 
notations in above theorem. In our case, the first task we must face is the estimate 
of the norm of U*'^~^{fj.)£{fi)U on some suitable contours. The estimate will be 
taken with respect to oo-norm, but this task is very tedious, we put the details in 
the appendix and give only the statement as below here. In the following lemma, 
the contour is chosen as: for any complex number fio , 

r„(^,) = {^:|^-/io|=0(^).} (13) 

in which we do not specify the radius since we shall obtain a uniform estimate. 

Lemma 3 Let jjlq be a zero of the transcendental equation 

fi sin /uvr — avr cos /ivr = 0, 

locating in a suitable neighborhood of some sufficiently large integer n, where an 
is a characteristic value of Qi. Let Tn{fJ.o) is chosen as ( \13i) with the same n. 
Then 

n 

Using Largrange inversion formula, it is not difficult to see that the zero of the 
transcendal equation 

/x sin fiir — air cos fin = 0, 

lying in a suitable neighborhood of positive integer n can be expressed as 

a K(a) 

fin = n-\ \ ^ H . 

n n"^ 

By Theorem |^ and Lemma |3|, if an is a characteristic value of Qi, then there 
exist at least one eigenvalue A„(q) of (||) whose square root vA„(a) satisfies 

\/%ia)=n + -+jnia), (14) 

n 

where 7n(a) = o(— j). We summarize above as 



Lemma 4 Suppose the eigenvalues of ^ satisfying ( |J^ in asymptotic sense. 
Then 



{air : all possible a that (I4) holds } 5 a{Qi), 
where o'{Qi) is the spectral set of Qi. 
On the other hand, we have the following result. 

Lemma 5 // any sequence of eigenvalue of ^) has the asymptotic expression as 
(jTjj, then a-K G cr{Qi) and 7n(a) = o(l/n^). 



Proof. By®, 

W(/i2) = -^sinnTTl + cos inrGi + ^^^^^^^02 + ^^^^^Gs + 0{-^). (15) 



2\ • r , r I sin^vr cos^vr 1 



Let fin = vA„(a) as given by (|T3 ), then plugging /i„ into (|T5[), we have 



W{l4i) = -{n-l |-7n(a))sin(nH h 7n(«))7r/ + cos(n H \- ^nM)'^Qi 

n n n 

1 / a / ^^ ^ 

-^ T^-^ ^^sin(nH h7n(a))7r^2 (16) 

ra + f + 7n(a) n 



- — — — —-2 cos(n + -7„(a))7rg3 + 0{— 



Since 



sin(//,7r) = (-l)-(^+^„(a))7r + 0(^), 



n n^' 



cos(/.„7r) = (-l)"(l-i(«+^„(a))V)+0(ij; 

2 n n^ 



Substituting them into (16) and equating the terms by comparing orders, we have 



1 



Wifii) = i-irianl -gi+ n7„(Q)7r/] + 0(— ). (17) 



Since W^fj.^) is singular for all n sufficiently large, and n7„(a) -^ 0, by the conti- 
nuity of determinant function, the matrix anl — Qi must be singular. Otherwise, 
letting T = avrl — ^1, we have, by (|l7|), that 

\\W{fil)-{-lTT\\<\n^n{a)A+0{\). 

Thus if T is invertible, then so is W{fi'^) for n sufficiently large, which is absurd. 
Therefore, avr is a characteristic value of Qi. 

By Lemma ^ and Lemma ^, we have Theorem ^ . 



Remark. Above approaches can be extended to some other type of boundary 
conditions, the author has disposed some kinds of them in |^], e.g., the Dirichlet 
boundary conditions. The determination on the first term in (^) is quite precise, 
while on the second or higher terms, there are still some ambiguities need to be 
excluded. 
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Appendix. The estimate of ||?7*^ ^{fi)S{fi)U\^_ 

In section 2, we give the statement on the estimate of ||C/*^^^(/i)i?(/i)C/||oo- The 
proof of that lemma lies on the estimate of orders on the prescribed contour. 
By (loD and iM), 



where 7?.(;u) = diag[^i(/Li), . . . , ^at (//)]. Let //o be a value making \I'(/io) be singular. 
Then fio is a zero of Qjin) for some j. According to the distribution of the 
zeroes of the transcendental equation Qj{ij,), we may assume fj,o lie in a suitable 
neighborhood of a sufficient large integer n, and can be espressed as 

a b ^/ 1 X 

/io = nH h — + 0{—), 

n Ti'^ n^ 

where 5 depends on a. With this kind ^o as center, we choose r„(/io) as defined 
in (|l^). We shall show that: On such contours, the two quantities 

,,^, .sinUTT,, ,,„, .COSUTT,, 

are of order 0(l/n), and hence the remainder term is of order less then 0(l/n). 

Here, we only give the demostration on the first quantity. As for the other 
one, one can use the arguements in the coming paragraphs to prove it. 

Write 

( xiip) ••• \ 

X2{^^) ■■■ '■■ 

: : ■•. 

V ••• XJv(/u) ) 

where Xi(/^) = sin fin / fj,{fi sin fin — ajircos fin), 1 < j < N. Then 

\\TZ{fi) lloo = max ||xi(Ai)||oo = max sup |xj(^)|. 

M i<J<^ i<j<Af^er„(Ma) 
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We only to give the estimate on |xj(/^)l ^^ the sense of growth of order. Write 
H = s + it, then for any /i G r„(/io), we have 

a Scoscj) dsincj) b ^/ 1 n 

fJ, = n-\ 1 ^ \-t ^ 1 ^ + C(— =■), 

n n^ n^ ri'^ nP 

where 5 is the radius of the contour and < </> < 2tt. Then, 

a 6 cos (f) b ^, 1 ^ , , 

n n^ m-^ n^ 

5 sin 6 , , 

t = -^. 19 

The fohowing identities are needed: 

|sin/x7r| = — (cosh 2t7r — cos 2s7r), 

I cos /xvr p = - (cosh 2t7r + cos 2s7r) , 

23f?{/i sin yuvr cos /ivr} = s sin 2s7r + t sinh 2t7r, 



and 



|//sin/x7r — ^vrcos^vrl = |/isin/i7r| 



+ (6*71)^1 cos/ivrp — 2^7r3f?{/isin/U7rcos/i7r}, 
where S R. Plugging (|l8|), (p!9|) into above identities, we have 



a^vr^ 27r2 



Isin/uvrP = — t, — I ^a(5cos</) 

7.2 1 

+ ^(-aV + 6a6 + 3<52)) + 0(— ), if a / , 

|sin/u7r|2 = ^^+ 54^_2(,os2(^) + o( )^ if a = 0. 

|2 1 «'^' 27r2 
cos^vr = 1 T. ^ao cos 

7.2 1 

+ ^(3<52(a^^2 _ g^5 ^ 1 _ 2cos2 </,)) + O(^), if a / 0, 

7^2 1 

Icos^uttP = 1 + — r(l-2cos2(/>) + 0(^), if a = 0. 



Hence we have : 

27r2 
|;usin^7r| + {Ott) |cos^7r| ={9 + a )tt -\ a5 cos (p 

Ji 1 

+ ^(3(5^ - aV + 6a^ - 'iijiaGf + 6a6) + O(^), if a / 0, 
3n^ n'^ 

l/isin^uvrp + (6'7r)2| cos^uvrp = (6'7r)^ H ^ + 0(^-), if a = 0. 

On the other hand, we have 

23f?{^ sin ;[/7r cos /uvr} = 2a7r H (5cos(/> 

n 

27r 1 

+ ^(3a2 + 36-2^2^3)^Q(_)^ ifa/0, 

27r 1 

23f?{;Usin;U7rcos/i7r} = — (5cos(^ + 0(— r), if a = 0. 

n n** 

Therefore, we obtained 

l/i sin yuvr — 07r cos [i'k\ 

27r2 7r2 

= (avr -d^f + b(a - 0) cos + ^ (6(a - Q) {c? + 6) 

n 3n^ 

+ 35^ - T?a^{a^ + 30^ - \aQ)) + O(^), if a / 0, 



l^usin/ivr — 07rcos;U7r| = (07r) 6'cos0H ^ — h 0(— r), if a = 0. 



|2 _ ^fl_^2 _ _ 
Finally, we have : On r„(/io). 



and 



sin^yrl = — |a|7r + o(— ), if a / 0, 
sin^vrl = — ^Jvr + of^r), if a = 0. 



l^sin^uvr — 07rcos/U7r| = 0(1), if 7^ a. 

I^sin^uvr — 07rcos/U7r| = — (Jtt + o( — ),if a 7^ 0, if = a. 

n n 

l/isinyuvr — 07rcos/U7r| = — (57r + o(— ), if a = 0, if ^ = a. 

n n 
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So, if a / 0, we obtain 



-lalvr + ofi) , 1 

(n + o(l))(|0-a| + o(l)) ~^n2 






|X,(^)I = 7 J«l^ + Q(n) =i^|c,vr|+o(l), ifg = a. 

' (n + o(l))(i^5 + o(i)) n^' 

while a = 0, 

^n^'' ~ 7 — ^ — l^\\l\a\ I — nvT ~ ~3"mT' ii ^ 7^ a, 
(n + o(l))(|t^| + o(l)) n-^ |t^| 

^57r + o(A) 1 1 

(n + o(l))(i(^7r + o(^)) n^ n^ 

Summarize above results as : for n sufficiently large, 

sin fiTT 1 I I 

||7c(/i) lloo ~ -^|a|7r, 1cm if a 7^ 0, 

/i on 

n{f^) ^ oo ~ 0K> ifa = 0. 

This tells us that ||'^(^) '^'°/'^ ||oo has growth order at most 0(l/n) as n is 
sufficiently large. With such estimate oon the contours, we have the desired 
estimate on WU*'^^^ {fi)£ {fj,)U\\, so we complete the proof Lemma ^ . 
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